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ABSTRACT 



Atmospheric density fluctuations induce phase 
perturbations that degrade the spatial coherence of a laser 
beam propagating through the atmosphere. These 
degradations spread the laser beam and alter the centroid 
and intensity profile stochastically. Turbulent conditions 
arising from various environmental situations and 
meteorological phenomena are found at virtually all levels 
of the atmosphere. A substantial fraction of the optical 
turbulence along a vertical path arises from the surface 
heat flux within the first 100-200 meters above the ground. 

This thesis seeks to measure and analyze these 
turbulent layers. A high frequency acoustic echosounder 
was developed to analyze atmospheric turbulence within the 
first 100-200 meters above the ground with extremely high 
accuracy. The echosounder design incorporated a 25 
element, square, planar array housed within a lead and foam 
insulated, acoustic shroud. This acoustic echosounder was 
employed to collect real-time, low level atmospheric data 
at two sites. This data should prove useful in determining 
whether the performance of laser and electro-optical 
systems could be increased by raising sensors and 
transmitters well above the local turbulent layer through 
the use of towers or other such structures. 
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i. i N mmuc.x i . ON 



A coherent laser beam propagating through the 
atmosphere is very susceptible to numerous turbulent- 
dependent processes. As electromagnetic waves transit the 
turbulent regions, atmospheric irregularities randomize the 
amplitude and phase of the wave. A high resolution 
acoustic profiler was designed and developed to quantify 
the altitude dependence of these atmospheric disturbances. 
Acoustic energy is more susceptible to these irregularities 
than electromagnetic waves. Historically acoustic sounders 
have proven to be the best means of providing high 
resolution atmospheric turbulence profiles within the first 
500 to 1000 meters above the ground. Acoustic echosounders 
are frequently used to detect and measure atmospheric 
density and velocity fluctuations resulting from wind 
shears, convection, and temperature inversions. 

Presently two optical parameters, the spatial coherence 
length ( ro ) and the isoplanatic angle (0o) are measures of 
the perturbation of electromagnetic waves propagating 
through the atmosphere, and are measured very reliably by 
optical systems developed by Walters [Refs. 1 and 2]. 
Although these systems measure a line integral of the 
atmospheric optical turbulence with high accuracy, a major 
drawback of each system is that no provision is provided 
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for determing the height of the atmospheric disruptions. 
If these disturbances are found to exist very near the 
surface, it may be possible to negate their effect by 
elevating sensors and transmitters or altering the surface 
heat flux to suppress the generation of turbulence. 

This thesis deals with the design, construction and 
implementation of a high frequency, acoustic echosounder 
which will accurately analyze the atmospheric density 
fluctuations within approximately the first 200 meters of 
the atmosphere. As this project is a product of the 
research and efforts of two students, the work was 
appropriately divided. My particular task was concentrated 
in the actual design and hardware development of the 
echosounder. My colleage, LT . Michael Wroblewski [Ref. 3], 
dedicated his efforts toward software development and data 
display, analysis and storage. 

The echosounder developed and discussed herein utilizes 
a 25 speaker, square, planar antenna array enclosed within 
an acoustically absorbent shroud comprised of lead and foam 
insulation inside a sturdy plywood framework. This antenna 
array was employed as both the transmitter and receiver of 
the acoustic signal. A computer controlled function 
generator transmits an acoustic pulse through the antenna 
array. Additional electronic components incorporated 
within the overall echosounder system are used to amplify, 
digitize, display, analyze and ultimately store the return 
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signal. Over time the received signal data is used to 
produce a quantitative, three dimensional map of the 
atmospheric turbulence. 

Acoustic echosounders have been developed and in use 
for many years and have proven to be a valuable probe for 
analyzing the structure and dynamics of the lower 
atmosphere [Refs. 4 through 7], Devices similar to ours 
have been used to obtain profiles of the atmospheric 
density and temperature fluctuations [Refs. 8 and 9], Our 
device excels by utilizing a high speed HP217 computer to 
control and monitor the echosounder. This enables us to 
obtain real-time data collection and analysis with the 
ability to store and reproduce the atmospheric profile 
plots at will. This information should help quantify the 
effects of low level turbulence on electro-optical systems 
performance. It may then be possible to negate these 
turbulent effects by simply mounting sensors and 
transmitters well above the turbulent layers or by adaptive 
optics . 

This thesis discusses the development of a high 
frequency, 5KHz , acoustic echosounder, the criteria 
involved in the actual designs of both the acoustic array 
and the acoustically absorbent enclosure, and the selection 
and function of all the associated component parts of the 
echosounder system. Finally the operational success of the 
system will be discussed and the output data analyzed. 
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II. BACKGROUND 



Acoustic echosounders probe the atmosphere by 
transmitting a pulse of acoustic power which is 

subsequently scattered back from the atmosphere by 

temperature and velocity inhomogeneities. The echosounder 
(echosonde) equation, often referred to as the radar 
equation in meteorology, is used to determine the 

backscattered acoustic power. This equation is summarized 
by Neff in Reference 10 and is based upon the work of 
Tatarski [Ref. 11] and Little [Ref. 5], 

Pr = Er[ PtEt ][ exp( -2aR) ][ ao(R,f) ][ ^ct ][ AGR- 2 ] 

where 

Pr is the electrical power returned from a range R, 

Pt is the electrical power transmitted at frequency f, 
Er is the efficiency of conversion from acoustic power 
to electrical power by the transducer, 

Et is the efficiency of conversion from electrical 
power to acoustic power by the transducer, 

exp(-2aR) is the round trip power loss due to 
attenuation where a is the average attenuation (meters - 1 ) 
to the scattering volume at the range R (meters), 

ao(R,f) is the scattering cross section per unit volume 
at a distance R and frequency f, 
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c is the local speed of sound (meters/second), 
t is the pulse length (seconds), 

A is the aperture area of the antenna (meters 2 ), 

R is the range (meters), and 

G is the effective-aperture factor of the antenna. 
Empirically measuring or calculating the values for all 
other terms, one can use this equation to determine 
0o(R,f), the scattering cross section per unit volume; that 
is, the fraction of incident power backscattered per unit 
distance into a unit solid angle at a frequency f. Based 
upon experimental results, Tatarski [Ref. 11] expresses the 
acoustic backscatter cross section per unit volume, 
cro ( R, f ) , by , 



k = 27r/A is the incident acoustic wavenumber at 

wavelength A , 

To is the local mean temperature in degrees Kelvin, and 
Ct 2 is the temperature structure parameter. 

Combining this equation with the echosonde equation, one 
obtains a volume-averaged measure of Ct 2 . 



Ct 2 



ao ( R , f ) = 0.0039 k 1 / 3 



To 2 



where 



1 



1 



2 



1 Pr 



Ct 2 



To 2 k-i/3 



R 2 exp ( 2aR ) 



0.0039 ErEt 



ct AG Pt 
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Hall and Wescott [Ref. 12] calculated a beam-shape 
compensation factor of 0.40 for a piston source antenna 
with a uniformly illuminated square aperture. This value 
is the same as the effective-aperture factor, G and can be 
substituted into the above equation. Approximating the 
aperture area of the antenna to be equivalent to 25 times 
the aperture area of a single speaker having a diameter of 
7.620 centimeters, we get a value of 0.1140 square meters. 
Combining these values with the numerical constants in the 
above equation, we can simplify the equation for Ct 2 . 



1 1 Pr 

Ct 2 = 11245 To 2 k-i/3 R 2 exp(2ccR) 

Er Et ct Pt 



The efficiency of conversion factors for our 



echosounder were determined with the use of an IVIE 
Electronics IE-30A Audio Analyzer Sound Level Meter, a 
HP3314A Function Generator, a HP3561A Dynamic Signal 
Analyzer, a reference speaker, the acoustic array and the 
Naval Postgraduate School anechoic chamber. The function 
generator was used to transmit a 1.0 volt, 5.0 KHz, 
sinusoidal signal through the reference speaker. The 
acoustic array was employed as a receiver and was aligned 
along the centroid of the main lobe of the reference 
speaker signal. With the acoustic array in position, the 
sound level meter was used to determine the total power 
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available for reception by the acoustic array; that is, the 
total power actually ensonified within the solid angle of 
the acoustic array. Ten independent measurements were made 
with the sound level meter and then averaged to yield a 
power value of 73.7 dB. Using a conversion factor and the 
array aperture area determined above, the total acoustical 
power, Pa, available at the acoustic array is given by the 
equation below. 

Pa = 107 . 37 . ( i . 10-12 Watts / m2 ) • ( 0.11401 m2 ) 

or, 

Pa = 2.67 • 10 -6 Watts 

The dynamic signal analyzer was then used to measure the 
actual power received by the acoustic array and converted 
into electrical power. Four independent measurements were 
recorded and averaged to yield a value of 4.035 mVrms. The 
acoustic array impedance was previously measured and 
determined to be 12.1 ±0.4 Q. These values were used in 
the following equation to determine Pe , the total 
electrical power converted by the acoustic array. 

( 4.035 • 10-3 Vr ms ) 2 

Pe = = 1.34 • 10- 6 Watts 

12. 1 Q 

The efficiency of conversion factor for acoustic to 
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electrical power is the ratio of Pa to Pe , or, 



1.34 • 10" 6 Watts 

Er = = 0.503 

2.67 • 10-6 Watts 

The value for the efficiency of conversion for electrical 
to acoustic power, Et , is not readily calculatable with the 
equipment at hand. However, based upon the design of the 
speakers used, the values for Er and Et should be 
approximately equal. Furthermore, allowing for 

experimental error especially in the reading of the sound 
level meter, the best accuracy for the values of Er and Et 
is 0.5 ±0.1. Substituting conversion factor values of 0.5 
into the above equation, we can further simplify the 
equation for Ct 2 . 

To 2 Pr 

Ct 2 = 44980 k~i/3 R 2 exp(2aR) 

CT PT 

The above equation was incorporated into a computer 
program that provides the system control, data acquisition 
and data processing techniques [Ref. 3], The reduced data 
was then used to measure the temperature structure 
parameter as a function of time of day and altitude for 
various sites. 
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III. SYSTEM DESIGN AND EQUIPMENT DEVELOPMENT 



A . SPEAKERS 

In the design of our echosounder, it was determined 
that a rapid decay time was required to obtain accurate 
short range information. Weight restrictions involved with 
equipment transportation plus high efficiency in the 2KHz 
to 5KHz range led to our decision to use piezo ceramic 
speakers. The Motorola KSN 1005A speaker was selected 
based on these requirements and the specifications charted 
in the Motorola catalog [Ref. 13] and reproduced below in 
Figures 1 through 3. 



Nominal Power, Impedance, and Distortion Ratings 




300 500 1000 2000 5000 10000 20000 40000 

Fig. 1. Speaker Ratings 
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Typical Frequency Response 




FREQUENCY (Hz) 

500 1000 2000 5000 10000 20000 40000 



Fig. 2. Speaker Frequency Response 



Dimensions: KSN 1005A, KSN 1003A 
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These charts indicate that a maximum response for our 
speakers occurs at a resonant frequency of 5000 Hertz. 
This frequency was used as the baseline from which all our 
measurements were made. 

Using the speakers in an anechoic chamber, the average 
e _ 1 voltage decay time was measured to be approximately 900 
usee (Fig. 4). This decay time translates into a sound 
propagation distance of just over 15.0 centimeters (at STP) 
from the speakers. Considering our requirements, this 
speaker was ideally suited to serve our purpose. 




Fig. 4. Speaker Decay Time Trace 
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B. ACOUSTIC ARRAY 



Our next. consideration was the echosounder array 
pattern. Ideally, the acoustic sources should be placed 
exactly one half wavelength apart. At a frequency of 5000 
Hertz, this would require spacings of 3.4 centimeters (at 
STP) which is physically impossible for the speakers we 
have chosen. The closest possible spacing was 7.62 
centimeters between sources after shaving off the flange of 
the horn (Fig. 3). 

From Kinsler, et al . [Ref. 14: p. 195], the equation 
for the directionality factor of a simple line array is 
derived as: 



H ( N , 9 ) 



1 

N 



N 

sin ( kd sin 9 ) 

2 



1 

sin ( kd sin 9 ) 

2 



where 

k is the wave number (2k/ \ ), 
d is the distance between sources, 

N is the number of sources, and 

9 is the angle measured from a line perpendicular to 
the array to the direction of interest. 

However, this equation assumes simple point sources 
which does not accurately describe the speakers we have 
chosen. It was necessary to couple this equation to the 
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directionality factor for a piston source which is also 
identified in Kinsler, et al . [Ref. 14: p. 108]. 



D(9) 



2 Ji ( ka sin 9 ) 



ka sin 8 



where 

k is the wave number, 

a is the radius of the piston source, 

8 is the angle measured from a line perpendicular to 
the array to the direction of interest, and 

Ji is a first order Bessel function. 

These two equations were combined to produce the 
equation of directionality for a linear array of piston 
sources, L(N,8), by simple multiplication. 

L( N , 8 ) = H( N , 8 ) • D( 8 ) 



This equation was incorporated into a computer program 
[Appendix A] and used to generate linear array patterns for 
varying numbers of piston sources (Figs. 5 through 9). 
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Fig. 5. Three Element Array Fig. 6. Four Element Array 




Fig. 7. Five Element Array 
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It was concluded from these plots of the linear array beam 
patterns for varying numbers of piston sources that five 
linear elements would provide the best combination of 
forward directionality, sidelobe suppression, physical 
size, and a relatively low cost. Then, in order to enhance 
both array efficiency and symmetry, we settled on a five by 
five element array design with 7.62 centimeter spacing 
between speakers in both the vertical and horizontal 
directions . 

After we verified the manufacturer’s polarity 
designation for 35 speakers and inspected each of these 
speakers for physical defects and loose wiring connections, 
we obtained Lissajous plots for each individual speaker. 
Figure 10 represents a typical Lissajous plot for one of 
the speakers incorporated into the five by five element 
array. Based upon the speaker’s output to input voltage 
ratios as illustrated by these Lissajous plots, we were 
able to rank all of our speakers by signal efficiency. It 
was based upon this criteria that we selected the 25 most 
efficient speakers for use in the five by five element 
array, placing the best speakers at the center and 
subsequent ranking speakers further toward the sides and 
corners . 
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Fig. 10. Lissajous Plot 

The 25 selected speakers were mounted in a five by five 
planar array on a balsa wood insulated bilayered sheet 
metal board. After wiring all the speakers in series, we 
surrounded all the electrical connections and speaker backs 
with two 3.0 centimeter layers of foam insulation 
sandwiching a 1.0 millimeter lead sheet. Then the entire 
array mounting was enclosed in a 44 by 44 by 5 centimeter 
sheet metal box. This design was chosen to suppress 
virtually all acoustic energy propagating out the rear 
hemisphere of the array, while shielding the array from any 
external electrical interference (Figs. 11 and 12). 
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Fig. 11. Array Photograph 
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C . ENCLOSURE 

Acoustic echosounding has proven to be an extremely 
useful technique for probing and analyzing the lower 
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atmosphere. In order to utilize the acoustic waves 
transmitted and later received by this remote sensing 
method efficiently, it is essential to have an efficient 
antenna with highly directive beams and strongly suppressed 
sidelobes. Antenna design becomes increasingly more 
important in a noisy environment where noise pollution 
within the sidelobes may dominate the desired signal within 
the main lobe. Hall and Wescott [Ref. 12] showed that 
sidelobe suppression improved with higher frequencies. 
Their studies showed that the measured 90 degree sidelobe 
suppression ranged from 38 dB at 1 KHz to 50 dB at 5 KHz. 
Furthermore, any significant improvement in sidelobe 
suppression could only be obtained by surrounding the 
antenna with an acoustic energy absorbing cuff or shroud. 

In an effort to maximize our antenna main lobe to 
sidelobe power ratio, we intend to operate at 5KHz . 
Additionally, we have designed an acoustic energy absorbing 
enclosure. Many designs were considered based upon 
previous research in the field of echosounding [Refs. 15 
through 18]. In addition, we obtained the actual acoustic 
beam patterns for our array using a computer program 
written by LCDR Butler [Ref. 19] which we modified for our 
purposes. This modified version of LCDR Butler’s computer 
program is in Appendix B. By rotating the array in an 
anechoic chamber, we were able to produce highly accurate 
polar plots of the array beam patterns (Figs. 13 and 14). 
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5X5 ACOUSTIC ARRAY BEAM PATTERN 




R=3. 60 M MIKE VT= 8.2465 INPUT’S. 0 V FREQ’S. 0 KHZ 



Fig. 13. Polar Plot of 5 X 5 Array Beam Pattern 
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DIAGONAL 5X5 ACOUSTIC ARRAY BEAM PATTERN 




Fig. 14. Polar Plot of 5 X 5 Array Beam Pattern 

at a 45° Aspect 



27 



A comparison of the theoretical beam pattern for the 
five element, piston source array in Figure 7 with the 
experimentally measured polar plot of the 5 X 5 acoustic 
array in Figure 13 and the associated coordinates 
corresponding to the individual data points illustrate 
significant agreement. As predicted by our computer model 
[Appendix A] , the actual acoustic beam pattern is comprised 
of nine distinct lobes in the forward hemisphere. Both the 
positions of these lobes and their relative amplitudes are 
within 2.0 percent of their predicted values. The computer 
modeled main lobe full width angle of 19.5 degrees very 
closely predicted the experimentally measured main lobe 
full width angle of 19.0 degrees. Reflected, scattered, 
diffracted and electronic noise all contribute to the 
clutter near the origin of Figure 13. This is the primary 
reason why the minima don’t actually go to zero as 
predicted by our model. Additional reasons include 
mechanical noise in the device used to rotate the acoustic 
array and vibrational noise in the array itself . The 
apparent lobe in the upper left hand quadrant is an anomaly 
seen in all polar plots herein, and corresponds to acoustic 
energy reflected off the door and window of the anechoic 
chamber. Figure 14 represents the acoustic array mounted 
at a 45° aspect. A computer model for this orientation is 
a very complicated task and was not performed at this time. 
However, this polar plot of the acoustic array beam 
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pattern, and in particular, the appreciable sidelobe 
suppression strongly support the need for further array 
design studies. 

Since it was our aim to suppress all sidelobes and 
utilize the main lobe, we chose not to taper our enclosure 
as most previous researchers had. Rather we designed the 
enclosure based upon the dimensions of the array itself and 
the acoustic beams it generated. The height of the 
enclosure is 160 cm. Allowing for the base foam and lead 
insulation, the array housing and the base framework 
itself, the actual array piston sources are positioned 143 
cm from the mouth of the enclosure. The length and width 
of the enclosure are both 65 cm. Again allowing for 
construction materials, the inner dimensions are both 48 
cm. Hence, the position of the center speaker is 24 cm 
from any of the four sides. These dimensions produce the 
unobstructed 19.0 degree path into which the main lobe is 
transmitted . 

2 • tan -1 ( 24 cm / 143 cm ) = 19.05® = 19.0® 

Plywood was used for the construction of the enclosure 
and provided not only a rigid, inexpensive framework, but 
also provided an impedance mismatch which helped to 
attenuate external noise interference. Anticipating all 
kinds of weather conditions during data collection, the 
plywood enclosure was first waterproofed with four coats of 
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marine varnish. Grooved joints, caulking and weather 
stripping were also design considerations. 

A millimeter layer of lead can suppress an acoustic 
signal by reflection as much as 40 dB (Figs. 15 and 16). 
About 7.0 centimeters of convoluted foam can suppress a 
signal another 3 to 4 dB (Figs. 15 and 17). The primary 



purpose of the 


convoluted 


foam is 


to 


attenuate the 


sidelobes which 


repeatedly 


reflect 


off 


the enclosure 



interior at glancing incident angles. Together the lead and 
convoluted foam make an extremely efficient absorbing 
material for use in our enclosure. 




Fig. 15. Signal Response Reference, No Insulation 
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START 4 000 Hz RANGE: -51 dBV STATUS: PAUSED 

A: MAG RMS : 25 




Fig. 16. Signal Suppression by Lead 



START 4 090 Hz 
A: MAG 



RANGE: -21 dBV STATUS: PAUSED 
RMS : 25 




X: 5000 Hz Y : -22 . 73 dBV 



Fig. 17. Signal Suppression by Foam 
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Lead-lined absorbing foam was glued to all inner 
surfaces of the enclosure with two 1.0 millimeter layers of 
lead overlapping at all corners. Strong aluminum brackets 
were used to connect the four side panels to each other as 



well as to the enclosure base (Figs. 18 through 20). 




Fig. 18. Enclosure Photograph, Fully Assembled 
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D. ADDITIONAL HARDWARE 



In addition to the acoustic array and enclosure 
discussed in the previous sections, there are primarily six 
more pieces of equipment essential to the function of the 
echosounder. Below is a brief description of these 

additional components . 

1 • HP200 Series Computer 

The Hewlett-Packard (HP) 200 Series Computer 
includes a 20 megabyte hard disk drive, a floppy disk drive 
and an associated printer and monitor. The HP200 Series 
Computer is the central control component for the entire 
echosounder arrangement. The computer used is an HP217 
programmed in Basic 3.0 and equipped with an Infotek BC203 
Basic Compiler and an Infotek FP210 Floating Point 
Accelerator to enhance the speed of program execution. The 
computer is used to maintain timing and process control, to 
define all parameters for the HP3314A Function Generator as 
well as to execute the trigger command that produces the 
echosounder transmitted acoustic pulse. The computer also 
processes the data from the acoustic array via the pre- 
amplifier, bandpass filter and analog to digital converter. 
The computer then conducts data reduction routines to 
produce and display a high resolution atmospheric profile. 

2 . HP3314A Function Generator 

The Hewlett-Packard (HP) 3314A Function Generator 
is a multimode function generator capable of providing 
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sine, square and triangular wave functions as well as any- 
desired waveform ranging in frequency from 0.001 Hertz to 
19.999 Megahertz. The HP3314A Function Generator was used 
to supply a pulse of an integer number of sinusoidal cycles 
of constant amplitude to the QSC Model 1700 Audio 
Amplifier. Typically a pulse of 100 sinusoidal cycles was 
used. A constant frequency setting of 5000 Hertz was used 
for all data runs. 

3 . QSC Model 1700 Audio Amplifier 

The QSC Model 1700 Audio Amplifier is a high power 
amplifier which can supply 350 watts over an 8 ohm load. 
This amplifier was used to boost the output voltage of the 
function generator by a factor of 20 from 1.5 volts to 30 
volts. The power supplied to each speaker in the array 
during operation was then 37.1 watts as indicated in the 
equation below. 

( Voltagerms )2 ( 30 )2 

Power = = 37.1 Watts 

Impedance 2 • 12.1 

4 . Pre-amolif ier 

The pre-amplifier was designed and constructed by 
Walters [Ref. 20] to supply a gain of 1000 to the returned 
signal and isolate the data aquisition system from the 
transmitted pulse. The pre-amplifier may be thought of as 
a safety and switching mechanism for the system. An LT 
1007 and an LT 1037 Operational Amplifier were selected for 
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use in the pre-amplifier based upon their low noise 
properties which were evaluated by physically incorporating 
them into the device and measuring their noise and gain 
characteristics with a HP3561A Dynamic Signal Analyzer. 
The IN 4000 Series Rectifiers isolate the power amplifier 
from the receiver between pulses and limit the voltage 
applied to the operational amplifiers. The electrical 
diagram of this device is outlined in Figure 21. 




Fig. 21. Schematic of Pre-amplifier 
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5 . Rockland Wavetek Model 852 Bandpass Filter 



The Rockland Wavetek Model 852 Filter operates as a 
48 dB per octave bandpass filter that suppresses broadband 
noise in the system. High and low bandpass settings of 
5500 Hertz and 4500 Hertz respectively were used with a 
roll off of 48 dB per octave. The filter response at these 
settings is illustrated in Figure 22. 



START 1 000 Hz RANGE s -5 dBV STATUS: PAUSED 

A: NAG RMS: 50 




START: 1 000 Hz BU: 477.43 Hz STOP: 31 000 Hz 

X: 5000 Hz Y:-20.24 dBV 



Fig. 22. Bandpass Filter Response 



The 3dB (half power) bandwidth was approximately 1 KHz 
wide. Although it would have been desirable to reduce the 
bandpass interval to between 50 and 100 Hertz without 
impairing signal response, the Rockland Wavetek Model 852 
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Bandpass Filter was the best filter available. It 
performed adequately despite the large 1 KHz bandpass. 

6 . Infotek AD200 Analog to Digital Converter 

An Infotek AD200 12 bit Analog to Digital Converter 
was used to sample the amplified signal voltage and 
digitize it for the computer system. A typical sampling 
frequency used was 20000 Hertz. 

The arrangement of the above equipment is illustrated ' 
in Figure 23. 




Fig. 23. Echosounder Equipment Arrangement 
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E . SOFTWARE 



The software and the HP217 Computer were responsible 
for controlling and monitoring every operational aspect of 
the system. LT Wroblewski [Ref. 3] provides an in-depth 
discussion of the echosounder software which is also 
summarized in Appendix C. 
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IV. DATA ANALYSIS 



A . ECHOSOUNDER PERFORMANCE 

An analysis of "the acoustic echosounder output was 
conducted to determine the validity of previously 
determined echosounder parameters. The e~ 1 decay time 
constant for the applied voltage of a single speaker was 
previously calculated to be 900 us which corresponds to a 
distance of 15.3 cm. The e~ 1 decay time constant for the 
applied voltage of the entire acoustic array was determined 
by the same experimental procedure to be 1060 us which 
corresponds to a distance of 18.0 cm. Although the 25 
speakers in the acoustic array transmitting in unison 
should have approximately the same e~ 1 decay time constant 
as the single speaker, it is not surprising that the array 
e - 1 decay time constant is larger. Using ten time 
intervals, an expectation of 99.99 percent decay of the 
transmitted signal should occur within a distance of 180.0 
cm or just under 2.0 meters. Reviewing typical echosounder 
traces and Ct 2 plots [Appendix D] , it is evident that such 
a short recovery time is actually quite accurate. The 
recovery distance is typically found to be on the order of 
9 meters or 53 ms. This "blind zone" is not uncommon in 
most radar systems and is in part a return signal from 
ground clutter surrounding the acoustic array and its 
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enclosure immediately following pulse transmission. An 
analog to digital intensity factor of 4000 is assigned to 
the heavily shaded region of the recovery zone. Since the 
computer software has generated this value as the product 
of the actual analog to digital "count" value and the range 
from the acoustic echosounder, this value of 4000 must be 
divided by the range at the recovery zone boundary to give 
us an actual analog to digital "count" value. This result 
is then multiplied by an experimentally measured conversion 
factor of 3.34 • 10 -8 voltsr ms /count . This result may then 
be used in the determination of the total electrical power 
received and converted by the acoustic array. The value 
for the total electrical power received by the acoustic 
array is calculated in the equations below. 



Voltager ms = (3.34 • 10 -8 v/count) ( 4000m» counts ) (9m) -1 

Voltagerms = 1.48 • 10 -5 volts 

Power = (Voltagerms) 2 • (Array Impedance) -1 
Power = (1.48 • 10 -8 volts) 2 • (12.1 Q) -1 
Power = 1.82 • 10 -11 Watts 



□sing this received electrical power value and the known 
transmitted acoustic power value of 37.1 Watts, we can 
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solve for "the e - 1 decay time constant of the applied 
voltage. The product of the total acoustic power received 
and the efficiency of conversion factor ( Er = 0.50 ), 
yields the total electrical power received. Hence, the 
total acoustic power received is 3.64 • 10" 11 Watts. 
Although the distance of 9 meters corresponds to a time 
interval of 53 ms, this value is actually inflated by the 
system pulse length of 20 ms. Compensating for this, the 
actual recovery zone time interval is 33 ms. The e - 2 decay 
time constant of the power is equal to the e - 1 decay time 
constant of the applied voltage. These values are 
calculated below. 



3.64 • 10-ii Watts 



37 . 1 Watts 



e - 2t /x or, 



x 



- 2 • 33 ms 
In (6.37 • IO -12 ) 



2.56 ms 



) 



This e ~ l decay time constant of the applied voltage 
corresponds to a distance of 43.5 cm and is about twice the 
previously determined value. 

Consideration was also given to the construction of a 
hexagonal planar array rather than a square planar array 
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based upon the astounding sidelobe suppression for the five 
by five element, square array at a 45° orientation as 
illustrated in Figure 14. To investigate the potential of 
such a close packed design, two small seven element 
hexagonal arrays were constructed. The first array design 
utilized a close packed spacing of 7.62 cm, while the 
second array design employed a 3/2 wavelength or 10.2 cm 
spacing. Polar plots of the beam patterns for both these 
hexagonal array designs were made by rotating the arrays in 
an anechoic chamber with the main diagonal both horizontal 
and vertical. In both cases the first sidelobe suppression 
for both close packed orientations was down by nearly 15 
dB, while the first sidelobe suppression for both 3/2 
wavelength spacing orientations was only down by about 7 
dB. Furthermore, the ratios of the size of the main lobe 
to those of the sidelobes is much greater for the close 
packed design than they are for the 3/2 wavelength spacing 
arrangement for both the vertical and horizontal 
orientations. The polar plots of the array beam patterns 
for both close packed and 3/2 wavelength spacing hexagonal 
designs are illustrated in Figures 24 through 27. These 
figures strongly warrant further study of the hexagonal 
close packed array design. 
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HORIZONTAL DIAGONAL, CLOSE PACKED, 7 SPEAKER HEXAGONAL ARRAY 




MIKE 



R-3. 60 M MIKE VT- 5. 6364 INPUT-5. 0 V FREQ-5. 0 KHZ 



Fig. 24. Polar Plot of Close Packed Hexagonal Array 

(Horizontal Diagonal) 
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VERTICAL DIAGONAL. CLOSE PACKED. 7 SPEAKER HEXAGONAL ARRAY 




MIKE 



R=3. 60 M HIKE VT- 5. 1559 INPUT-5. 0 V FREQ=5. 0 KHZ 



Fig. 25. Polar Plot of Close Packed Hexagonal Array 

(Vertical Diagonal) 



46 



HORIZONTAL DIAGONAL. 3/2 WAVELENGTH SPACING. 7 SPEAKER HEXAGONAL ARRAY 




MIKE 



R-3.60 M MIKE VT- 5. 18 INPUT-5.0 V FREQ-5.0 KHZ 



Fig. 26. Polar Plot of 3/2 Wavelength Hexagonal Array 

(Horizontal Diagonal) 



47 



VERTICAL DIAGONAL. 3/2 WAVELENGTH SPACING. 7 SPEAKER HEXAGONAL ARRAY 




Fig. 27. Polar Plot of 3/2 Wavelength Hexagonal Array 

(Vertical Diagonal) 
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B. SITE EVALUATION 



Echosounder data was collected at two different 
locations. The primary data collection site was the upper 
roof of Spanagel Hall at the Naval Postgraduate School, 
Monterey, California. This site was chosen simply for 
convenience of data collection. Data gathered at this 
location is believed to represent the California coast 
during the spring near sea level. The second data 
collection site chosen was in the vicinity of the 24 inch 
telescope at Lick Observatory, San Jose, California. This 
site is located atop Mt. Hamilton at an altitude of 
approximately 5700 feet and nearly 20 miles inland from the 
coast . 

These two data collection sites represent areas of 
differing atmospheric air pressures, water vapor pressures, 
local temperature ranges, and local wind velocity ranges. 
These characteristics all play important roles in effecting 
the local atmospheric turbulent conditions and thereby the 
atmospheric structure parameter, Ct 2 . 

In addition to collecting echosounder data at Lick 
Observatory, simultaneous measurements of the isoplanatic 
angle (9o) and spatial coherence length ( ro ) were made with 
systems developed by Walters [Refs. 1 and 2]. A basic 
knowledge of these two systems is necessary to understand 
the correlation procedures made. The isoplanatic angle 
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( 9 o) is primarily an upper atmospheric measurement which 
indicates atmospheric disruptions at a range of two to 15 
kilometers above the Earth’s surface. On the otherhand, 
the spatial coherence length (ro) is a uniform measure of 
the effects of the entire atmospheric blanket on coherent 
light transmission. A close comparison of the acoustic 
echosounder data, the isoplanatic angle ( 0 o ) measurements 
and the spatial coherence length ( ro ) measurements should 
give us an indication of the turbulent contributions from 
the lower and upper troposphere as well as the stratosphere 
above . 

Based upon isoplanatic angle ( 9o ) and spatial coherence 
length (ro) measurements at Mt . Wilson near Los Angeles, 
California, a strong correlation between the isoplanatic 
angle ( 80 ) and spatial coherence length ( ro ) measurements 
occurs if the low altitude boundary layer contribution is 
sufficiently small. This strong correlation helps to 
reinforce the overall description of the atmosphere at the 
time of data collection. Figures 28 and 29 graphically 
illustrate the atmospheric measurements made at Mt. Wilson 
on 2 April 1987 by Walters. The strong correlation between 
the isoplanatic angle (9o) and the spatial coherence length 
( ro ) is especially evident between the hours of 0700 and 
1300 universal time. The close tracking of these two 
measurements during this time interval indicate that the 
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Fig. 28. Isoplanatic Angle Measurements, Mt. Wilson 




Fig. 29. Spatial Coherence Length Measurements, Mt. Wilson 
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upper atmospheric conditions, as measured by the 
isoplanatic angle ( 9o ) , are dominating the entire 
atmospheric profile as measured by the spatial coherence 
length fro). Unfortunately, the existence or non-existence 
of any turbulent surface effects cannot be ascertained by 
the employment of the above two systems alone. However, 
employment of these two systems together with the acoustic 
echosounder should enable us to produce a complete 
atmospheric profile with strong correlation between all 
three atmospheric measurements. 

On 9 and 10 April 1987, all three systems were operated 
at Lick Observatory. Again good correlation between the 
isoplanatic angle ( 9o ) and the spatial coherence length 
( ro ) measurements was noted. In addition, strong 
correlation between the spatial coherence length ( ro ) 
measurements and the acoustic echosounder data was 
conclusively present. A comparison of the isoplanatic 
angle ( Qo ) and the spatial coherence length ( ro ) 
measurements in Figures 30 and 31 shows good correlation 
between the two atmospheric parameters especially during 
the 0930 to 1130 time interval on 9 April 1987. However, 
during subsequent hours the isoplanatic angle ( Uo ) 
measurements remain high (approximately 12 u.rad ) indicating 
relatively calm turbulent conditions in the upper 
atmosphere, while the spatial coherence length ( ro ) values 
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drop sharply after 1100 universal time indicating dominant 
increasing lower atmospheric turbulence. This increase 
in the low level turbulence should be evident in the 
echosounder data commencing around 1200 universal time 
(0400 local standard time) on 9 April 1987. A comparison 
of the echosounder data in Figures 32 through 34 
illustrates this increase in the local surface turbulence. 




Fig. 30. Isoplanatic Angle Measurements, Lick Observatory 
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Fig. 31. Spatial Coherence Length Measurements, 

Lick Observatory 
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Fig. 32. Echosounder Trace, Lick Observatory 
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Fig. 33. Echosounder Trace, Lick Observatory 



9 Rpr 1987 LICK OBSERVRTORY 




Fig. 34. Echosounder Trace, Lick Observatory 
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The strong correlation between the echosounder data around 
0400 to 0500 local standard time and the spatial coherence 
length ( ro ) measurements at 1200 universal time, combined 
with the lack of correlation between the isoplanatic angle 
(00) and the spatial coherence length ( ro ) measurements 
indicate that the lower atmospheric and surface turbulent 
layers are dominating the atmospheric profile during this 
time period. 

Data collected on 10 April 1987, again illustrate the 
strong correlation between the three atmospheric 
measurements made. A comparison of the isoplanatic angle 
(0o) and the spatial coherence length ( ro ) measurements 
during the time interval of 0700 and 1300 universal time 
indicate steady turbulent conditions in the upper 
atmosphere and greatly varying turbulent conditions at 
lower atmospheric levels. This is evident in Figures 35 
and 36 by the consistent values of the isoplanatic angle 
(0o) during the time period compared with the steady 
increase and eventual decline of the spatial coherence 
length ( ro ) values during the same time interval. The 
trace variations in Figure 36 during the hours of 0800 and 
1200 universal time are indicative of a period of 
diminishing lower atmospheric or surface turbulence 
followed by the onset of an increasingly turbulent period 
around 1200. This turbulent trend is distinctively 
supported by the echosounder data in Figures 37 through 44. 
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Fig. 35. Isoplanatic Angle Measurements, Lick Observatory 




Fig. 36. Spatial Coherence Length Measurements, 

Lick Observatory 
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RRNGE IN METERS RRNGE IN METERS 



10 Rpr 1387 LICK OBSERVATORY 




Fig. 37. Echosounder Trace, Lick Observatory 



10 Rpr 1387 LICK OBSERVATORY 




Fig. 38. Echosounder Trace, Lick Observatory 
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Fig. 39. Echosounder Trace, Lick Observatory- 
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Fig. 40. Echosounder Trace, Lick Observatory 
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Fig. 41. Echosounder Trace, Lick Observatory- 
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Fig. 42. Echosounder Trace, Lick Observatory 
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Fig. 43. Echosounder Trace, Lick Observatory 
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Fig. 44. Echosounder Trace, Lick Observatory 
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The large 200 mm spatial coherence lengths around 1100 
universal time are consistent with the low turbulence 
evident in the echosounder profiles around 0300 local 
standard time. 

Only acoustic echosounder measurements were made on the 
upper roof of Spanagel Hall at the Naval Postgraduate 
School. Data runs on 26 and 27 April 1987, include both 
echosounder data and the associated Ct 2 plots. Of 
particular interest is the plot of 11:30 on the 27th of 
April in Figure 45, which shows the maritime boundary 
inversion layer around 100 m being perturbed by convective 
plumes at lower altitudes. Sixteen additional echosounder 
output traces were included in Appendix D to exhibit 
typical atmospheric activity. Many of the plots, such as 
the 14:45, 15:30, and the 17:30 of 26 April and the 10:00, 
11:15. 13:00, 14:15, 15:15, 16:15, and the 17:15 of 27 
Apri L , have convective plumes which are prevalent whenever 
a heat flux between the surface and atmosphere exists. The 
plots of 18:15 and 19:15 on the 26th of April and of 18:00, 
18:45, and 20:00 on the 27th of April clearly show the 
passage of the neutral event which is encountered when the 
atmospheric and surface temperature difference becomes 
negligible. Finally, the plot of 16:30 on 26 April can be 
associated with strong winds which exhibit a somewhat 
uniform return for all altitudes across the entire 15 
minute interval of the trace. 
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V . CONCLU SIONS AND RECOMMENDATIO NS 



A high frequency, 5kHz, phased array echosounder 
constructed to measure low level turbulence appears to work 
within the 10 to 135 meter altitude range. When using the 
device, detailed profiles of the short range atmospheric 
density fluctuations were obtained. The profiles were 
found to correlate very well with the measurements of the 
isoplanatic angle (9o) and the spatial coherence length 
( ro ) during periods of simultaneous operation. This short 
range echosounder, when used in conjunction with the other 
atmospheric measuring devices, is an invaluable tool . The 
graphical output provides us with a more complete 
description of the atmosphere that can be used to calculate 
Cn 2 , the atmospheric index of refraction structure 
parameter . 

Further software developments such as the incorporation 
of Doppler and Fast Fourier Transform routines into the 
echosounder program will determine the radial velocity 
profile of the return signal. Such an improvement will 
determine the velocity of the probed air masses and plot 
these echosounder traces as a function of color intensity. 

Other areas of further reasearch include the design and 
testing of different array patterns. One such pattern, the 
hexagonal array, is already undergoing tests. As the 
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speaker array evolves, it is also evident that the 
enclosure must follow suit to accomodate the new array 
design and associated beam patterns. Further improvements 
in software routines are also inevitable. It is always 
desirable to store on a disc all data collected at a site. 
Presently, the floppy disc capacity is inadequate for the 
storage of more than two hours of data. Furthermore, if 
FFT and Doppler routines are to be added, the computational 
speed of the computer will inhibit the pulse repetition 
rate. A special data acquisition technique called "Direct 
Memory Access" may then have to be introduced into the HP 
computer to allow simultaneous data collection and 



processing . 

Finally, this product may be used in various 
applications not already addressed. One such use may be to 
measure windshear at airports. Other researchers have 
expressed interest in using the echosounder to measure 
arctic atmospheric conditions during meteorological 
surveys. Assuredly, as the device evolves into a more 
compact and highly efficient instrument, its range of 
application will continue to expand. 
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APPENDIX A 



LINEAR ARRAY PROGRAM 

! ANTENNA BEAM PATTERN FOR AN N ELEMENT LINEAR ARRAY 

GCLEAR 

Ri= . 038 1 

Fr=5000 

Ka=( Ri*2*PI*Fr )/343 

INPUT "ENTER THE # OF DRIUERS n ,N 

D= . 0762 

Kd=(2*PI*D*Fr)/343 

Kd2-Kd/2. 

Nurc_decades=4 
Scale=1 0 A Nuro_decades 
Thetastep=.001 
GCLEAR 

UIEWPORT 20,80,20,30 
WINDOW -20,105,-2,1 10 . 

AXES 1 00/ ( 2*Num_decades ) ,100/(2. *Nim_de cades ) ,0 ,0 ,20 ,20 
MOUE 10,103 

LABEL "PATTERN FOR ";Ns“ ELEMENTS" 

CSIZE 3 
FOR 1=1 TO 5 

MOVE -1 5 , ( 25* ( 1-1 )>-1 
M=( 10-( (6-1 )*10 ) ) 

LABEL M 
NEXT I 
LOIR PI/2 
MOUE -15,40 
LA8EL "DECIBELS" 

LOIR 0 

FOR Theta=Thetastep TO PI/2 STEP Thetastep 
Arg=Ka*SIN( Theta) 

Arg2=2*FNJ 1 ( Arg > 

Arg2=Arg2/Arg 
Kds=Kd2*SIN( Theta ) 

Htheta=SIN(Kds*N)/(N*SIN(Kds ) ) 

Htheta=Htheta*Arg2 
Htheta=Scale*Htheta*Htheta 
Htheta=LGT(Htheta )/Num_de cades 
X=Ht he ta*SIN( Theta ) 

Y=Htheta*COS( Theta ) 

PLOT X* 1 00 ,Y* 1 00 
IPRINT Theta, Htheta 
NEXT Theta 
END 
! 
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450 DEF FNJKX) 

4G0 ! J1 BESSEL FUNCTION 10 MAY ,1984; DLW 

470 IF X>3 . THEN GOTO Xgt3 
480 X3=X/3 . 

490 X3=X3*X3 

500 Fx = . 5+X3*< - . S6249985+X3*( . 2 1 093S73+X3*( - . 039S4289+X3* 

( .0044331 9+X3*( -.00031761 +X3*. 00001 109 ) ) ) ) ) 

510 J 1 =Fx*X 
520 RETURN J1 
530 Xgt3: X3=3./X 

540 Fx=.79788456+X3*( . 00000 1 56+X3* ( . 01 6S9667+X3*< . 000 1 71 05+X3* 

(-.0024951 1+X3*( .001 136S3-X3*. 00020033 ) ) ) ) ) 

550 Tx=X-2 . 356 1 9449+X3* ( . 1 2499G 1 2+X3*( . 00005G5+X3*( - . 00G37879+ 
X3*( . 00074348+X3*( . 00079824-X3* . 00029 1 G6 ) ) ) ) ) 

560 J 1 =Fx*COS( Tx )/SQR( X ) 

570 RETURN J1 
580 FNEND 
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APPENDIX B 



LCDR BUTLER’S BEAM PATTERN PROGRAM 



10 PRINT "** BEAM PATTERN COMPUTATION AND PLOT PROGRAM »*" 
20 PRINT "*■***■*■*■*■*■***■*■****■***■*■*******■**■***■**■*■*****■***■****" 
30 ! 



40 ! 

50 PRINT "***** SET_UP ***•*" 
60 ! 

70 ! 



80 PRINT CHRS (10) 

90 PRINT "THIS PROGRAM CONTROLS THE FOLLOWING APPARATUS" 
100 PRINT "TO MEASURE. DISPLAY. PLOT AND STORE THE" 

110 PRINT "DIRECTIONAL CHARACTERISTICS OF A SINGLET," 

120 PRINT "DOUBLET. TRIPLET, AND QUADRUPLET ARRAY." 

130 ! 



140 


PRINT 


" INSTRUMENT 




ADDRESS " 


150 


PRINT 


M 




• i 


160 


PRINT 


" HP 3478 MULTIMETER 


(VOLTS) 


724 


170 


PRINT 


" HP 3478 MULTIMETER 


(RES) 


723 


180 

190 


PRINT 

| 


" DATA AQUISIT ION/CONTROL 


709 



200 DISP "INITIALIZE INSTRUMENTS" 
210 ! 



220 REM *** PRESET AND INITIALIZE INSTRUMENTS** 
230 OUTPUT 709 ;"RS" 

240 OUTPUT 723 
250 OUTPUT 724 



! RESET 

"F3R5T 1 " ! FREE RUN 30K0HMS 
"F2RAT1" ! FREE RUN AUTORANGE AC VOLTS 



260 ! 

270 DISP "INITIALIZE PROGRAM" 

280 ! 

290 OPTION BASE 1 

300 COM ANGLE ( 150) ,CALC_ANGLE( 150) , NAMES [80 3 
310 DEG 



320 PRINTER IS 701 



330 INTEGER I.C.L.N 

340 DIM R( 360 ) .VOLT (360) ,OHMS_DATA(360) ,X(360) ,Y<360) ,K(36G) 
350 MVOLT-O 

360 MASS STORAGE IS " . DATA1 " 

370 DISP "NUMBER OF DATA POINTS DESIRED?" 

380 BEEP 10,300 
390 INPUT N 

400 DISP "ENTER DISC DATA FILE NAME" 

410 INPUT NAMES 

420 DISP "WANT NEW DATA (Y/N)?" 
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430 BEEP 10.300 
440 INPUT 0$ 

450 IF Q$="Y" THEN GOTO 480 
460 IF Q$="N" THEN GOTO 920 
470 GOTO 430 

480 REM *** DEFAULT PARAMETERS *** 

490 ! 

500 ! MUST INSERT FOR LINES 390, 400, 410 AND 810 
510 ! PRIOR TO START OF RUN 
520 ! 

530 CREATE NAMES, 360, 16 
540 ASSIGN* 1 TO NAMES 
550 NAMES="RUN_1 " 

560 FULLSCALE=500 ! LOCK_IN MVAC 

570 CONV=360/1 0000 ! CONVERT OHMS TO DEGREES 

580 SET_UP : ! *** SET_UP MEASUREMENT PARAMETERS *** 

590 PRINT CHRS (10) 

600 PRINT ”1. GRAPH TITLE : M ;NAMES 

610 REM *** MODIFY IN PRECEEDIMG LINE WHEN NECESSARY *** 
620 DISP "MANUALLY ENTER DESIRED FREQ FOR SPEAKER DRIVER" 
630 DISP "ENTER 0 IF NO CHANGE" 

640 DISP "ENTER 1 TO CHANGE TITLE" 

650 BEEP 10,600 
660 INPUT C ! CHOICE 
670 IF C=0 THEN GOTO 690 
680 IF C=1 THEN GOTO 880 
690 PRINT CHRS (10) 

700 OUTPUT 709 ;"CLS0" ! ROTATE CW 

710 FOR 1=1 TO N 

720 OUTPUT 709 ;"CLS1" ! POWER TO MOTOR ON 

730 WAIT .5 ! MILLISECONDS 
740 OUTPUT 709 ;"OPN1" ! POWER TO MOTOR OFF 
750 WAIT 4000 

760 ENTER 723 ; OHMS_DATA < I ) 

770 ENTER 724 ; VOLT(I) 

780 PRINT "I IS NOW"; I 

790 PRINT "OHMS_DATA=" ;GHMS_DATA< I) 

800 PRINT "VGLT = " ; VOLT < I ) 

810 NEXT I 
820 FOR 1=1 TO M 

830 PRINT* 1 : OHMS_DATA( I) , VOLT ( I) 

840 NEXT I 

850 ASSIGN* 1 TO * 

860 GOSUB PLOT_ROUTINE . 

870 GOTO 910 

880 ENTER NAMES ! NEW TITLES 
890 GOTO 600 
900 PRINT "DONE" 
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S10 END 

920 DISP "WHAT IS DATA RUN DESIRED?" 

930 BEEP 10,500 
940 INPUT NAMES 
950 ASSIGN# 1 TO NAMES 
960 FOR 1=1 TO N 

970 READ# 1 ; OHMS_DATA( I ) , VOLT < I ) 

980 NEXT I 

990 ASSIGN# 1 TO * 

1000 GOTO 860 
1010 ! 

1020 ! 

1030 PLOT_ROUTINE : ! *** PLOTS POLAR DATA *** 

1040 ! 

1050 ! 

1060 PLOTTER IS 705 
1070 GRAPH 
1080 GCLEAR 
1090 PEN 1 

1100 LIMIT 50,210,20,180 

1110 SCALE -60,60,-60,60 

1120 XAXIS 0,10,-40,40 

1130 YAXIS 0,10,-40,40 

1140 CSIZE 3 

1150 LINE TYPE 1 

1160 FOR L=0 TO 40 STEP 10 

1170 MOVE 0 ,L 

1180 LORG 6 

1190 LABEL L-40 

1200 NEXT L 

1210 LINE TYPE 1 

1220 CSIZE 5 

1230 MOVE 0,50 

1240 ! 

1250 ! INSERT LABELS 1250 AND 1280 PRIOR TO START OF RUN 
1260 ! 

1270 LABEL ;"5 X 5 ACOUSTIC ARRAY BEAM PATTERN" 

1280 CSIZE 3 
1290 MOVE 0.-45 

1300 LABEL ;"R=3.60 M'V'MIKE VT = " , VOLT ( I ) , "INPUT = 2.5 V"." 
" , "FREQ = 5 . 0 KHZ" 

1310 CSIZE 4 
1320 MOVE 45,0 
1330 LORG 5 
1340 LABEL ;"MIKE" 

1350 LINE TYPE 4 

1360 FOR L= 1 0 TO 40 STEP 10 
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1370 MOVE L.O 

1330 FOR 1=0 TO 360 STEP 5 
1390 DRAW L*COS <I>,L*SIN <I> 

1400 NEXT I 
1410 NEXT L 
1420 LINE TYPE 1 
1430 PEN UP 
1440 PEN 2 

1450 FOR 1=1 TO N ! BEGINS BEAM PATTERN COMPUTATION 
1460 IF VOLT < I ) >MVOLT THEN MVOLT = VOLT ( I) 

1470 NEXT I 

1480 FOR 1=1 TO N 

1490 K=OHMS_DATA( I) / 10000*360 

1500 K 1 = K 

1510 IF OHMS_DATA< 1X10 THEN K = K1+360/295 

1520 IF OHMS_DATA(I)>1 0500 THEN GOTO 1610 

1530 VECTOR=20*LGT <VOLT< I) /MVOLT) +40 

1540 LINE TYPE 1 

1550 IF VECTOR< 1 THEN VECTOR=0 

1560 PRINT I, K. VECTOR 

1570 X(I)=VECTOR*COS <K) 

1580 Y < I) =VECTOR*SIN <K) 

1590 PLOT X< I) , Y< I) 

1610 NEXT I 
1615 PEN UP 
1620 RETURN 
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APPENDIX C 



ECHOSOUNDER SOFTWARE 

The driving force of this system is found within the 
HP200 Series Computer. The software and computer system 
were responsible for controlling and monitoring every 
operational phase of the hardware components. The program 
which accomplishes this task was entitled "ACRDR" and is 
listed in Reference 3. This program was modeled after a 
program written by Walters [Ref. 20], but each program 
performs a distinctly different computational task. This 
program was written in HP Basic 4.0 and was compiled by an 
Infotek BC203 Basic Compiler to enhance the speed of 
execution . 

The program "ACRDR" is easily broken down into a number 
of blocks and subroutines which perform specific 
operations. These sections are outlined in a flowchart 
(Fig. 46). Each block is straightforward in its purpose, 
and the program is designed to be as helpful to the user as 
possible. As a prologue to the actual program code, there 
is a listing of all the program variables with a short 
desription of their use. Such a listing familiarises a 
user with the computations to be made and also provides 
quick references for any future modifications to be 
implemented . 
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The initialization procedures which set all the 
parameters necessary for data collection follow the 
variable definition section. These parameters control such 
things as the contrast between background and return signal 
and the setting of the computer’ s internal clock. All 
options of the HP3314A Function Generator are decided by 
the user at this phase. Another initialization involves 
preparing the Infotek AD200 Analog to Digital Converter for 
operation. Finally, internal arrays are dimensioned, 
integer variables are defined and the computer function 
keys are redefined to suit the "ACRDR" program. 




Fig. 46. Flowchart of Fundamental Program Operations 
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Now that the system is ready to operate, certain 



constants must be determined based upon atmospheric 
conditions at the time of data collection. From the 
ambient temperature, a speed of sound in air is determined 
by the equation in Kinsler, et al , [Ref. 14: p. 106]. 

C = Co [ 1 + (T/273) ]i/2 



where 

C is the local speed of sound, 

Co is the speed of sound in air at 0° Celcius, (331.6 
m/ s ) , and 

T is the local temperature in degrees Celsius. 

Using this value and a user decided sampling frequency and 
data point total, the maximum range of the traces is 
determined by the relation, 

RANGE = C/2 • l/( SAMPLE FREQUENCY) • N , 



where 

RANGE is the range of the echosounder in meters, 

SAMPLE FREQUENCY is the user desired rate of data 
acquisition, and 

N is the number of data points per trace. 

Typically, the values of many of the parameters are not 
varied. The value of N above is almost always taken to be 
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the allowed maximum of 16301 data points and the speed of 
sound can be safely estimated to be 340 m/s at room 
temperature. By defining a typical sampling frequency of 
20000 Hz , a range of approximately 135 meters is obtained. 

Another parameter determined is the atmospheric 
attenuation coefficient. This calculation is done in a 
subroutine obtained from Reference 21 and is necessary in 
this program for the calculation of the temperature 
structure parameter, Ct 2 . Other computations made prior to 
program execution include the pulse length, wavenumber, and 
the D-C offset of the equipment. 

At this point, the computer is finally ready for data 
collection. The screen setup displays a distance versus 
time plot and the internal clock of the computer is 
synchronized with the pulse from the HP3314A Function 
Generator. Following transmission of the pulse, the return 
signal is received, digitized, stored in an array and data 
reduction commences. A block averaging technique is used 
in which a block of data points the size of the number of 
emitted cycles is summed and averaged. Returns are plotted 
as darkened areas on the echosounder traces with the 
intensity of the darkened area directly related to the 
magnitude of the return signal. This same procedure is 
used in the computation of CT 2 with some very important 
differences. Aside from simply correcting for irr 2 
spherical divergence, Ct 2 is corrected for electronic gain. 
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the ratio of power returned to power transmitted, the 
efficiency of the speakers as transmitters and receivers, 
the area of the speaker array, the atmospheric attenuation, 
the pulse length, the temperature, the scattering cross 
section per unit volume at a specific range and frequency 
and finally the effective aperture factor of the antenna. 
Additionally, Ct 2 (Range) is averaged for a particular 
altitude over 15 minute intervals. Finally, after each 
pulse is reduced, a corresponding mean square noise level 
is determined. This is done by averaging each block 
average at maximum range until at least ten values have 
been used in the average. An upper limit of five over the 
average noise level (this corresponds to voltage 
fluctuations on the order of 10" 7 volts) is set on the 
rountine to avoid averaging any strong return signals ox~ 
anomalies such as passing aircraft. After the ten values 
are averaged, every subsequent pulse is averaged into all 
the preceding noise levels and removed from each subsequent 
return signal. 

Ultimately, after each 15 minute interval, or at the 
user’s request, the display terminal image is printed. At 
this time the Ct 2 computations are conducted and plotted on 
the screen and printer as range versus logic of Ct 2 . Upon 
completion of the printing of these plots, the system 
begins the data processing for the next 15 minute interval. 



76 



There are certain options built into the program to 
allow a user to change various aspects of operation. The 
function keys allow the user to change the local 
temperature, sample frequency or intensity factor during 
program execution. Additionally, the user can quit or 
restart the program, print the partial trace on the screen 
or elect to save a trace on a floppy disc. The save 
routine is invoked for the subsequent 15 minute interval 
after the appropriate function key is depressed. Saving a 
future trace may be ackward, especially if a user would 
like to keep an interesting trace which is presently on the 
terminal. This problem cannot be readily solved unless 
each trace is recorded to disc without user intervention. 
At present, this is not done because only eight traces ( two 
hours of data) can be written to a floppy disc before it is 
full. 
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